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A  generalised  solution  for  s  class  of  printed  circuit  antennas,  excited  by  s 
strip  transmission  line  is  presented.  The  strip  transmission  line  may  be 
embedded  laslde  or  printed  on  the  substrate.  As  an  example,  microstrip  dipoles 
elactromagnatically  coupled  (paras iticelly  excited)  to  embedded  strip  transmiasio 
line  have  been  analysed  accurately  and  design  graphs  are  provided  for  a  specific* 
substrata  material.  These  graphs  permit  the  establishment  of  e  design  procedure 
which  yields  tha  microstrip  dipole  length,  overlap,  offset,  and  substrate  thick¬ 
ness  with  the  goal  of  a  desired  input  match. for  a  given  substrate  material.  The 


ON  THE  MODELING  OF  ELECTROMAGNETI CALLY  COUPLED 
MICROSTRIP  ANTENNAS  -  THE  PRINTED  STRIP  DIPOLE+ 


BY 

P.B.  Katehl  and  N.  G.  Alexopoulos 
Electrical  Engineering  Departnent 
University  of  California  Los  Angeles 
Los  Angeles,  California  90024 


^Performed  under  U.S.  Army  Research  Contract  DAAG  29-83-K-0067  and  U.S.  Navy 
Contract  N00014-79-C0856  Mod.  P00005. 


mi 

3SI 


Mi 


ABSTRACT 


A  generalized  solution  for  e  class  of  printed  circuit  antennas,  excited  by 
a  strip  transmission  line  is  presented.  The  strip  transmission  line  may  be 
embedded  inside  or  printed  on  the  substrate.  As  an  example,  microstrip  dipoles 
electromagnet ically  coupled  (parasltically  excited)  to  embedded  strip  transmission 
line  have  been  analyzed  accurately  and  design  graphs  are  provided  for  a  specific 
substrate  material.  These  graphs  permit  the  establishment  of  a  design  procedure 
which  yields  the  microstrip  dipole  length,  overlap,  offset,  and  substrate  thick¬ 
ness  with  the  goal  of  a  desired  input  match  for  a  given  substrate  material.  The 
method  accounts  for  conductor  thickness  and  for  arbitrary  substrate  parameters. 
Comparison  with  experiment  shows  excellent  agreement. 


I  INTRODUCTION 

Thi*  paper  addresses  the  problem  of  rectangularly  shaped  and  electronag- 
netieally  coupled  printed  circuit  antennas.  In  this  case,  the  excitation 
mechanism  la  provided  by  a  strip  transulsslon  line  eabedded  Inside  the  substrate 
(see  Figures  1,  2)  which  couples  energy  parasltically  to  the  microstrip  antenna. 
Eleetroaagnet leal ly  coupled  microstrip  dipoles  [l]-[5]  have  been  investigated 
by  empirical  or  approximate  analysis  techniques,  and  a  very  approximate  model 
has  been  derived  for  the  rectangular  microstrip  patch  (6],  Since  the  radiation 
mechanism  of  a  microstrip  dipole  is  very  similar  to  that  of  a  microstrip 
patch  [2],  the  model  developed  In  this  article  is  applicable  to  the  analysis  and 
design  of  microstrip  elements  which  are  rectangularly  shaped  but  with  a  width 
■■•ller  than  the  element  length.  Reference  to  Figures  1  and  2  assists  in  emphasiz¬ 
ing  that  the  parameters  of  the  problem  are  arbitrary  in  the  development  of  the 
model.  Including  substrate  thickness  and  relative  permittivity  to  account  for 
previously  recognized  substrate  effects  (7]>[15).  In  addition,  the  thickness  of 
metallic  conductors  is  included,  strip  transmission  line  and  microstrip  antenna 
widths  may  differ  and  the  effects  of  the  microstrip  antenna  overlap  and  offset 
with  respect  to  the  transmission  line  on  the  current  distribution  is  investigated. 
Finally,  the  model  presented  in  this  paper  can  resolve  (see  Figure  3),  by  the 
proper  choice  of  parameters,  the  problems  of  the  electromagnetlcally  coupled 
rectangular  microstrip  patch  antenna,  the  electromagnetlcally  coupled  microstrip 
dipole,  the  microstrip  fed  rectangular  patch,  radiation  by  an  open  circuited 
microstrip  transmission  line,  as  well  as  variations  of  these  geometries. 

The  emphasis  of  this  article  is  centered  on  the  application  of  the 
theoretical  model  to  electromagnetlcally  coupled  microstrip  dipoles,  comparison 
with  published  experimental  results  [5]  and  provision  of  design  graphs  for  a 
specif ic  substrate  material.  The  model  is  developed  by  obtaining  the  current 
distribution  along  the  dipole  and  transmission  line,  including  sll  mutual 


Interactions.  The  Method  of  Moments  is  applied  for  the  determination  of  the 

V 

current  distribution  in  the  longitudinal  dimension,  while  the  current  dependence 
in  the  transverse  direction  is  chosen  so  as  to  satisfy  the  edge  condition  at 
the  effective  width  location.  Upon  determining  the  current  distribution,  trans¬ 
mission  line  theory  is  Invoked  to  evaluate  at  the  chosen  reference  plane,  the 
self  impedance  Zg  of  the  microstrip  dipole.  This  leads  to  a  design  procedure 
which  yields  the  microstrip  dipole  length,  overlap  and  offset  so  that  a  desired 
input  match  can  be  achieved  for  a  given  substrate. 

II.  ANALYTICAL  FORMULATION 

A.  Current  Distribution  Evaluation 


The  microstrip  dipole  and  strip  transmission  line  considered  here  have 
the  cross  sections  shown  in  Figure  2.  It  is  assumed  that  t^  and  t2  ere  very 
•mall  compared  to  wavelength  in  the  dielectric  (*1**2  <<  *®  *1®° 

assumed  that  the  currents  along  the  dipole  and  the  feeding  line  flow  on  the 
bottom  surfaces  Sd  and  S{.  The  electric  field,  under  the  assumption  made,  is 
given  by  Pocklington's  Integral  Equation 


E1^ 


E.  Jt 


G^(r/'r>)  *  Jv(?')ds’  (i  -  1,2) 


where  E*(ri  is  the  electrric  field  in  medium  i  (i  ■  1  air,  i  ■  2  dielectric)  and 
cj<?/?')  is  the  Dyadic  Green's  Function  in  medium  i  due  to  source  v  (v  •  1  strip 
dipole,  v  •  2  transmission  line).  The  Dyadic  Green's  function  is  given  by  the 


expression 


cj<r/?’)  -  [kjT  •  pj(?/r*) 
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Since  the  widths  of  the  dipole  end  alcrostrip  line  ere  fractions  of  the  wavelength 
in  the  dielectric,  it  can  be  aseuaed  that  the  currents  are  unidirectional  and 
parallel  to  the  x  axis.  Therefore,  the  current  vector  in  Equation  (1)  can  be 
written  la  the  fora 


Jv(^)  »  X  J^x'^y’) 

where  J*(x')  *n  unknown  function  of  x'  while  j^(y')  is  of  the  fora 
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Here,  we  is  the  effective  strip  width  given  by  w^  »  wV  +  26V.  The  paraaeter 

i  la  the  excess  half  width  and  it  accounts  for  fringing  effects  due  to  conductor 

thickness.  Interpretation  of  the  choice  for  the  current  density  dependence  in 

y'.  Indicates  that  the  edge  condition  is  satisfied  at  y'  -  ±wV/2,  which  is  an 

v  e 

equivalent  strip  of  rero  thickness.  At  y’  -  ,  the  current  density  remains 

finite,  as  is  the  case  for  nonzero  thickness  conductors.  This  choice  of 
dependence  in  y'  for  the  current  density  yields,  as  it  will  be  shown  very  accurate 
results  for  alcrostrip  dipole  resonant  length.  Formulae  for  effective  width 
exist  in  the  literature  [16], 117]  and  they  have  been  adopted  in  this  model  formu¬ 
lation.  The  concept  of  effective  width  is  thrust  upon  this  formulation,  since 
the  current  dependence  in  y'  is  assumed  and  not  derived  by  numerical  solution  of 
Pockllngton's  Integral  Equation  in  both  the  x  and  y  directions.  Equations  (1) 


and  (2)  are  combined  now  to  yield  (ds*  ■  dx'dy') 
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with  the  dyadic  F^(?/^)  given  by 


F^Cr/?1)  -  (?/?•)  n  +  F^x(^)« 


Substitution  of  (8)  Into  (7)  yields  the  following  expression  for  the  x  component 


of  the  electric  field: 


Bt-C  pffipty 


where 
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with 
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Further  definitions  necessary  for  Che  Interpretation  of  Equations  (9)  and  (11) 


are  provided  as: 
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■  Kronecker's  delta 


*  Kronecker's  delta 


In  addition, 

.  ,  .  ..2  .2.1/2 

bs  bs "  t2  *  Uo  “  fX  ‘  kl3 

fj(X,h)  ■  uosinh(uh)  +  ucosh(uh) 
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and  p  -  [(x-x')2  +  (y  -  y')2]1^2 
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The  zeros  of  f^A.h)  and  f2(A,h)  lead  to  TE  and  TM  surface  wave  modes 
respectively  [7]- [15]. 

In  order  to  solve  Equation  (8)  for  the  currents  and  J2,  the  method  of 
moments  is  used.  The  strip  dipole  and  fee  line  are  divided  into  +  1  segments 

with  v  -  1,2  respectively.  The  x'  dependence  of  the  current  is  written  as 

Nv 

J^(x’)  -  £l^f(x')  (15) 

x  n»l  n  n 


where  the  expansion  functions  f^x')  have  been  chosen  to  be  piecewise  sinusoidal 


functions  given  by 


V»*> 
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n-J. 
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sink! 


X  ,  <  x'  <  X 
n-1  —  —  n 


X  <  x'  <  X 
n  —  —  n+1 


(16) 


0  anywhere  else 

with  the  length  of  each  subsection.  By  substituting  Equation  (15)  into  (8), 


the  electric  field  can  be  written  as 


\^1,2  n-1 


-w  /2 


(17) 


Furthermore,  the  electric  field  is  projected  along  the  x  axi'-  (y  -  y*t  z  -  z1, 
i  *  1,2)  using  as  weighting  functions  the  basis  functions  (Galerkln's  Method) 
and  therefore.  Equation  (17)  reduces  to  a  matrix  equation  of  the  form 


O  [1^1  «  [V*] 
mn  n  m 


i  -  1,2 


1,2 


(18] 


where  [V*j  is  the  excitation  vector  of  order  N,  +  K,,  [1V]  is  the  vector  of 
m  x  *■  n 


unknown  coefficients  of  order  N,  +  N_  and  [Z  )  is  the  Impedance  matrix  of  order 

1  *■  mn 

(Nj  +  Nj)  x  (N^  +  Nj)  whose  elements  are  given  by 

dy* 
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(i9: 


Using  this  form  for  the  evaluation  of  the  elements  of  the  impedance  matrix,  one 
can  solve  the  matrix  equation  shown  in  (8)  to  find  the  unknown  coefficients  for 
the  current. 

B.  Self-Impedance  Evaluation 

For  the  evaluation  of  the  mlcrostrlp  dipole  self  impedance,  transmission 
line  theory  is  applied  and,  because  of  this,  unimodal  behavior  of  the  field  far 
from  the  atrlp/dlpole  coupling  region  is  essential.  In  order  to  satisfy  this 
requirement,  the  transmission  line  is  kept  very  close  to  the  ground  plane  giving 
a  ratio  wldth/llne-to-ground  distance  >  2.0. 


For  the  particular  geometry  considered  here,  a  unlmodal  field  is  excited 
under  the  transmission  line  and  as  a  result  the  current  distribution  beyond  an 
appropriate  reference  plane  forms  standing  waves  of  a  TEM-like  mode.  For  this 
reason,  the  mlcrostrlp  line  is  approximated  by  an  ideal  transmission  line  of 
characteristic  Impedance  Zq  which  is  terminated  to  an  unknown  Impedance  Zg 
(see  Figure  4) .  This  Quasi-TEM  mode  has  a  wavenumber  6  and  a  standing  wave 
ratio  SWR  equal  to  the  average  values  evaluated  from  the  original  current. 

If  the  origin  of  x  coordinate  is  taken  at  the  position  of  Zg,  then  the 
voltage  and  current  waves  on  the  ideal  transmission  line,  with  respect  to  this 
plane  of  reference,  are  given  by: 


V(x)  -  Ae“'3Bx  +  BeJBx 


and 


I(x) 


i-  [A.-*6x  -  B.16*] 
o 


Considering  the  positions  x  and  x  .  of  two  consecutive  maxima  and  minima, 

uhx  mm 

Equation  (20a)  gives 


I  1  -  ?  e  ■ax 

max  m  A _ 

B  326*min 
'.in  1  *  A  e 


Since  the  absolute  difference  between  x  and  x  is  equal  to  one-fourth  of 

max  min 

the  wavelength  of  the  guided  wave  "  xBin  "  ^g^^»  then  Equation  (21) 


can  be  written  as 


_  J2Bx 
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The  reflection  coefficient  T  is  defined  as 


and  by  considering 


1  *j? 

max  J2 
-  e 
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min 


then  Equation  (22)  gives  the  following  expression  for  the  reflection  coefficient. 


T(o) 


.I^ie 

Y  +  1 


-j2Bx 
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From  (20a)  and  (20b)  an  expression  for  the  self  impedance  Zg  evaluated  at  the 
position  xq  can  be  written  as  follows: 


Z  (x  )  1  +  T(x  ) 

s  o  _ q__  . 

z  i  -  r(x  )  K 

O  O 

If  now  this  plane  of  reference  is  considered  at  the  position  of  the  first 
current  maximum  d  from  the  end  of  the  open  circuited  transmission  line,  then 
the  self  impedance  measured  with  respect  to  this  plane  is  given  by 


Z  (d  ) 
s  max 


1  +  T(d  ) 
max 

1  -  T(d  ) 
max 


.  1  +  r(0)e 


J2Bd 


max 


1  -  T(0)e 


j2Bd 


max 


(27] 


Since  for  the  evaluation  of  T(d  ),  method  of  moments  together  with  transmission 

max 

line  theory  has  been  applied.  Equation  (27)  Indicates  that  the  evaluated  self- 
impedance  depends  on  the  characteristics  of  the  substrate  (cr,h)  as  well  as  on 
the  embedding  distance  of  the  transmission  line,  the  overlap,  the  offset  and 
the  length  of  the  dipole. 

The  fundamental  design  procedure  is  now  revealed.  One  wants  to  choose, 
for  a  given  substrate,  the  right  position  of  the  dipole  so  that  optimum  resonance 
or  in  other  words  perfect  match  is  obtained.  The  condition  of  optimum  resonance 
is  characterized  by  the  relation 


which  combined  with  (27)  gives 


T(d  )  -  0 
max 


T(x)  -  0 


This  means  that  one  wishes  to  find  the  geometry  which  perfectly  matches  the 
dipole  to  the  transmission  line. 

C.  Numerical  Results 

The  self  impedance  is  evaluated  using  Equation  (27)  as  a  function  of 
the  length  of  the  dipole  and  the  overlap  < 


ovp  S. 


where  is  the  part  of  the  surface  of  the  dipole  which  is  over  the  transmission 

line  (Fig.  4).  The  real  and  imaginary  parts  of  the  normalized  self  Impedance 

are  plotted  (Fig.  5)  for  a  substrate  thickness  h  ■  0.077",  dielectric  constant 

er  »  2.53  and  embedding  distance  bs  ■  0.0485".  The  width  of  the  Btrips  is 

w  -  0.060"  and  the  thickness  t  ■  0.00025".  The  length  of  the  dipole  is  varying 

between  0.373"  and  0.349"  and  the  overlap  takes  values  between  38%  and  86%. 

Figure  5  implies  that  there  exists a  particular  overlap  for  which  the  curve  of 

Z 

the  dipole  lengths  goes  through  the  point  ^  ■  1  or  in  other  words,  for  the 

o 

geometry  considered  above  there  is  one  value  for  the  overlap  and  a  specific 
dipole  length  so  that  the  dipole  is  perfectly  matched  to  the  transmission  line. 

It  is  of  Interest  to  investigate  now  the  dependence  of  the  normalized  self- 
impedance  on  the  offset  and  the  dipole  length.  For  the  geometry  used  previously, 
the  overlap  is  kept  at  50%  while  the  dipole  length  is  varying  between  0.373"  and 
0.349"  «nd  the  offset  takes  the  values  0*00",  0.02",  0.04"  and  0.08".  From 
Figure  6,  one  can  see  that  there  is  a  specific  offset  for  which  the  curve  of 
the  dipole  lengths  crosses  the  axis  of  the  normalized  self  resistance  at  the 
point  1.  Therefore,  for  this  value  of  the  offset,  there  is  a  dipole  length 
which  can  give  perfect  match. 


However,  as  It  was  mentioned  before,  the  normalised  self  Impedance  Is 
also  a  function  of  the  substrate  characteristics  and  the  position  of  the  trans-r 
mission  line.  In  order  to  examine  this  dependence,  the  dielectric  constant, 
width,  thickness,  overlap  and  the  distance  of  the  transmission  line  from  the 
ground  plane  are  kept  constant  to  the  values  above,  while  the  dipole  length  Is 
varied  again  between  0.373"  and  0,349",  and  the  embedding  distance  takes  the 
values  0,036",  0.048",  0,054"  and  0,066".  Figure  7  Indicates  that  there  Is 
a  unique  value  for  the  embedding  distance  which  combined  with  the  appropriate 
dipole  length  can  give  perfect  match. 

It  Is  also  very  interesting  to  see  how  the  current  distribution  changes 
as  a  function  of  the  embedding  distance  b#  around  optimum  resonance.  Figure  8 
displays  the  current  distribution  normalised  to  the  Incident  current  on  the 
transmission  line  for  h  -  0.099"  (und er coup led ) ,  h  *  0.093"  (perfectly  matched 
and  h  •  0.082"  (overcoupled).  The  substrate  Is  durold  (er  -  2.53),  the  dis¬ 
tance  of  transmission  line  to  the  ground  plane  Is  0.0285"  and  the  length  of 
the  dipole  1^  »  0.353".  As  it  is  shown  In  the  figure  the  current  on  the  dipole 
takes  Its  maximum  value  when  It  is  perfectly  matched. 

D.  Comparison  with  Experimental  Results 

The  theoretical  analysis  of  a  printed  strip  dipole  electromagnetlcally 
coupled  to  an  embedded  microstrip  line  is  tested  by  comparing  theoretical  results 
to  experimental  ones.  Stern  and  Elliott  [5]  measured  experimentally  the  self 
Impedance  of  strip  dipoles  with  rounded  comers  (Fig.  9)  printed  on  durold 
boards  (er  -  2.35)  with  substrate  thickness  h  -  0,077"  and  excited  by  a  micro- 
strip  transmission  line  In  the  dielectric  at  a  distance  from  the  ground  plane 
equal  to  0,0285".  The  width  of  the  strips  Is  w  •  0.06"  and  the  thickness 
t  ■  0.00025".  The  self  Impedance  was  measured  for  different  dipole  lengths 
and  its  normalized  values  are  shown  on  a  Smith's  chart  (triangles)  In  Fig.  9. 


The  solid  line  corresponds  to  the  theoretical  result.  If  one  defines  as  resonant 
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length  Lr  the  length  of  the  dipole  which  gives  s  pure  real  self  impedance, 
then  the  experimental  resonant  length  is  about  0.390"  while  the  theoretical 
one  0.379".  Therefore,  there  is  a  difference  2.751  with  approximately  2% 
resulting  from  the  different  shape  of  the  dipoles  which  were  studied  analytically 
(rectangular  atrlps)  and  measured  experimentally  (strips  with  round  comers)  [18] 
The  difference  between  the  theoretical  and  experimental  values  of  the  self 
impedance  is  due  to  two  reasons:  1)  different  shspe  of  the  dipoles,  11)  the 
fact  that  in  the  theoretical  evaluation  of  the  current,  the  hybrid  nature  of 
the  modes  propsgating  in  the  microstrip  wss  taken  into  account  while  for  the 
experiments  only  an  equivalent  TEM  mode  was  measured. 

III.  CONCLUSIONS 

An  effective  method  has  been  presented  for  the  accurate  modeling  of 
rectangularly  shaped  printed  circuit  antennas  which  are  excited  parasitically 
by  a  strip  transmission  line  embedded  in  the  substrate.  The  various  parameters 
Involved,  such  s  microstrip  antenna  and  transmission  line  dimensions  and 
location  as  well  as  substrate  permittivity  and  thickness,  are  completely 
arbitrary  In  the  model.  The  current  distribution  is  computed  by  the  method 
of  moments  in  the  longitudinal  direction  of  the  antenna  and  the  transmission 
line,  while  in  the  transverse  direction  it  is  chosen  as  the  Maxwell  distri¬ 
bution,  thus  satisfying  the  edge  condition  at  the  printed  element  effective* 
width.  Comparison  with  experiment  indicates  that  this  model  provides  excellent 


accuracy. 


FIGURES 


Figure  1: 

Side  and  top  view  of  a  printed  atrip  dipole 
excited  by  a  transmission  line  embedded  in 
the  dielectric. 

Figure  2: 

Cross  section  of  a  printed  atrip  dipole 
excited  by  a  transmission  line  embedded  in 
the  dielectric. 

Figure  3: 

Some  mlcrostrlp  antenna  geometries  resolvable 
by  proposed  model. 

Figure  4: 

Current  amplitude  on  the  atrip  dipole  and 
microstrip  transmission  line. 

Figure  5: 

Z  /Z  as  a  function  of  ic  and  L. . 

s'  o  ovp  J. 

Figure  6: 

Z  /Z  as  a  function  of  b  and  L. . 

ao  ax 

Figure  7: 

Z  /Z  as  a  function  of  the  offset  and  L  . 
so  1 

Figure  8: 

Current  amplitude  on  the  atrip  dipole  and 
transmission  line. 

Figure  9: 


Comparison  of  theoretical  to  experimental 
results. 
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Figure  2 

Cross  section  of  e  printed  strip  dipole  excited 
by  a  transmission  line  embedded  in  the  dielectric 


Current  amplitude  on  the  atrip  dipole  and 
microstrip  tranamiaeion  line. 
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